Local entanglement and quantum phase transition in the Hubbard model 
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The local entanglement E v of the one-dimensional Hubbard model is studied on the basis of its 
Bethe-ansatz solution. The relationship between the local entanglement and the on-site Coulomb 
interaction U is obtained. Our results show that E v is an even analytic function of U at half-filling 
and it reaches a maximum at the critical point U — 0. The variation of the local entanglement 
with the filling factor shows that the ground state with maximal symmetry possesses maximal 
entanglement. The magnetic field makes the local entanglement to decrease and approach to zero at 
saturated magnetization. The on-site Coulomb interaction always suppresses the local entanglement. 
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Quantum entanglement, as one of the most intrigu- 
ing feature of quantum theory, has been a subject of 
much study in recent years, mostly because its nonlocal 
connotation^ is regarded as a valuable resource in quan- 
tum communication and information processing 0, 0- 
For the application purpose, much of recent attentions 
have been focused on entanglement relevant to realistic 
systems. For example, several authors have investigated 
entanglement in spin systems 0, |M IE 0J1j El as weu as 
indistinguishable-particle systems [HJ, [HI • The work of 
Osterloh et o/.0 and Osborne and Nielsen^ on the XY 
model suggestively showing that the entanglement of two 
neighboring sites displays a sharp peak either near or at 
the critical point where quantum phase transition under- 
goes. Investigating the critical entanglement between a 
block of continuous spins and their supplemental parts in 
a spin chain model, Vidal et al. |9( pointed out its rela- 
tion to the entropy in conformal field theories. Recently, 
we studied entanglement and quantum phase transition 
of the XXZ model ^3 an d obtained its dependences on 
the anisotropy parameter A and correlation length £. 

It is well known that the Hubbard model is a typical 
model describing correlated Fermion systems. It plays 
a crucial role for understanding many physical phenom- 
ena in condensed matter physics, such as magnetic or- 
dering, Mott-insulator transition, and superconductiv- 
ity, etc. Therefore, the study of the entanglement in the 
Hubbard model not only provides possible clues for ex- 
perimental realization, but also sheds new light on the 
understanding of quantum many-body systems. 

In this Letter, we study the local entanglement of one- 
dimensional Hubbard model0. We obtain the depen- 
dence of the local entanglement on Coulomb interaction, 
particle number and magnetization respectively. We find 
that the ground state with maximal symmetry possesses 
the maximal local entanglement, and the Coulomb in- 
teraction always suppresses the local entanglement. To 
the best of our knowledge, no one has investigated the 
entanglement in the interacting many- fermion systems. 
Our results, which are based on the exact solution of the 
Hubbard model, will be inspirable for people to explore 
quantum entanglement and phase transition via nonper- 
turbative approach for other interacting many-fermion 



systems. An important observation of our studies is that 
the local entanglement is a smoothly continuous function 
at the quantum phase transition (Mott-transition) point, 
which clarifies that the singularity of the entanglement 
at quantum phase transition point is not an universal 
behavior as suggested by previous studies. 

The Hamiltonian of the one-dimensional Hubbard 
model reads 



n j] n ]i (1) 



, L, c - and c 1(T create and anni- 



where a =t, I; j = 1, ■ ■ . , L, c) a 

hilate respectively an electron of spin a at site j. There 
are four possible local states, |0)j, | T)j, I I T-L)j de- 
noted by \f)j, v = 1,2,3,4. The Hilbert space of L- 
site system is of 4 L dimensional, and \v\, v-z ■■■vr) = 

11^=1 \ v i)i are its natural basis vectors. Therefore any 
state in such a system can be expressed as a superposi- 
tion of the above basis vectors. We consider the ground 
state |\E') of the model. The local density matrix is a re- 
duced density matrix pj = TVj| 1 E r )(\E r |, where Tr., stands 
for tracing over all sites except the jth site. This en- 
ables us to study the correlations between the two parts 
of the Hilbert subspace defined on C 4 (g> (C 4 ) L_1 . Ac- 
cordingly, the von Neumann entropy E v calculated from 
the reduced density matrix pj is employed to measure 
the entanglement of states on the jth site with that 
on the remaining L — 1 sites. It was called the local 
entanglement for it exhibits the correlations between 
a local state and the other part of the system. 

The Hamiltonian possesses U(l)xSU(2) symme- 
try, i.e., it is invariant under the gauge transformation 
Cj.er - > e Cj.ct and spin rotation Cj. a — > U a sCj,s, which 
manifest the charge and spin conservation. Thus, so- 
lutions of Q should simultaneously be the eigenstates 
of particle number N = X)j( n i T + n j,i) an d the z- 
component of total spin S z = yV(re.,-,j — n j,i)- This 
implies the absence of coherent superposition between 
configurations with different eigenvalues of rij f and 
7ij.f — n^j.. The density matrix pj therefore has the fol- 
lowing diagonal form 

Pj = z\ o)(o | +u+ ItXTl +u~ ||)(|| +w IT1XT1I (2) 



Parameters in the density matrix pj represent local pop- 
ulations of four states, i.e., 

w = (n j]nji ) = Tr(n 3] -n nP: j), 
u + = (n-f) — w, u~ = {ni) — w, 
z = 1 — u + — u~ — w = 1 — (rij) — (nj.) + w (3) 

where (n-f) and (nj_) are electron density with spin up 
and spin down respectively. Thus the key point in dis- 
cussing the local entanglement is to study the density 
distribution of the four local states. 

As the Hamiltonian is invariant under translation, the 
local density matrix pj is site independent. Consequently, 
the von Neumann entropy (or the local entanglement 
which we call hereafter) of this local reduced density ma- 
trix is 

E v = —z log 2 z — u + log 2 u + — u~ log 2 u~ — w log 2 w 
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FIG. 1: The density of double occupancy w (left figure) and 
the local entanglement E v (right figure) versus the on-site 
coupling U . The solid lines are obtained from Eq. (|SJ and 
Q, while the open dots are obtained by solving the Eqs. (QJ 
of N — 70 sites system numerically. 



The one-dimensional Hubbard model has been solved 
exactly by the Bethe-ansatz method Q, 0] . The energy 
spectra is given by E = — 2 J2j=i cos kj, for U > and 
N < L. The charge and spin rapidities {kj, X a } fulfill the 
following transcendental equations: 

M 

2irlj = kjL — 9i(\ a — sinfcj), 

a=l 

N M 

2nJ a = 0i(\ a - sin kj) - ^ 6 2 {\ a -X b ), (4) 

j=l b=l 

where n (k) = 2 tan _1 (4fc/n<7), M is the number of elec- 
trons with down spins, and J a play the role of quantum 
numbers. The ground state is a spin singlet, whose solu- 
tion consists of real fcs and As, which is solved by choosing 
the quantum number configuration as successive integers 
or half-odd-integers symmetrically arranged around zero. 
For cases of either U < or N > L the energy spectra 
can be obtained by a particle-hole transformation 14J 

E(N - M, M; U) 

= (N- M)U + E(N — M,L — M; -U), 

= -(L-N)U + E(L-N + M,L-M;U) (5) 

Now we are in the position to investigate the rela- 
tionship between the local entanglement and the on-site 
Coulomb coupling U, the chemical potential, and the ex- 
ternal magnetic field. 

Local entanglement at half-filling and quantum phase 
transition driven by the on-site Coulomb coupling: ^From 
the Licb-Wu solution of the ground state at half-filling 
|14|. we obtain the density of double occupancy w by 
taking the derivative of Eq(U)/L with respect to U ac- 
cording to the Hellman-Feynman theorem, 

[°° M^J^duj 

W J a 1 + cosh([/ W /2) • ( ) 



where Jo(^) and J\(ijj) are zeroth and first order Bessel 
functions. The ground state is a spin singlet which im- 
plies (n^) = (?i|) = 1/2, hence u + and u~ are both 
simply 1/2 — w. The local entanglement is then 

E v = -2w log 2 w - 2 (1/2 - w) log 2 (1/2 - w) . (7) 

By making the use of Eq. (J3J), one easily finds that 
w(—U) = 1/2 — w(U), so the local entanglement is an 
even function of U 

E v (-U) = E V (U). (8) 

In the strong coupling limit, U — > oo, the double oc- 
cupancy tends to zero and the sites are singly occupied. 
The local density matrix pj has two nonvanishing values 
u + = u~ = 1/2 that results in E v (oo) = 1. For finite pos- 
itive U the hopping process | T)j| — * I ti)j|0)j+i 
or |0)j| tl)j+i brings about the occurrence of double oc- 
cupancies. The lower the on-site repulsive coupling U 
is, the larger the occurrence likelihood will be. Conse- 
quently, decreasing U will enhance the local entangle- 
ment E v . For negative U case, however, double occu- 
pancy is favored. If U — > — oo, half of the total sites 
are doubly occupied while the other half are empty, 
which gives rise to z — w — 1/2, u + = u~ = 0, and 
E v {— oo) = 1. For finite negative U, the hopping pro- 
cess | Tl)j|0)j+i -> | T)jl l>i+i or | l)j\ t) i+ i breaks up 
some of the double occupancies and suppresses the mag- 
nitude of w and z. This enhances the local entanglement, 
E v . As a result, the local entanglement increases when 
\U\ decreases. It is expected to reach a maximum value 
when U approaches to zero from both positive and neg- 
ative sides. We plot the density of double occupancies 
w and local entanglement E v as functions of U in Fig. 
H obtained by integrating Eq. JHJ. Data points in the 
figure are obtained by solving Bethe-ansatz equation (0J 
numerically on a 70-site chain. The excellent agreement 
justifies the validity of later calculations. 
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It is well known that the ground state of the one- 
dimensional Hubbard model at half-filling is metallic for 
U < 0, and insulating for U > 0. The point U = sep- 
arates metallic and insulating phases. Our calculation 
shows that the local entanglement reaches the maximum 
at the critical point U = where quantum phase tran- 
sition occurs. It is also worthwhile to notice that the 
derivative of the local entanglement with respect to the 
on-site coupling U is negative in the insulating phase but 
positive in the metallic phase. 

Eq. JJ]J could be integrated by performing series ex- 
pansion in U for the Bessel functions. In the large 
U ^> 1 region, to the third order in 1/U 2 , we have 
w = 4 In 2/U 2 - 27C(3)/[/ 4 + 375C(5)/[/ 6 where C stands 
for the Riemann zeta function. This shows that w de- 
scends rapidly when U increases. Therefore the local 
entanglement yields the following asymptotic behavior 



1 + 16\nU/U 2 



(9) 



Whereas in the week coupling region < U -C 1, 
the density of double occupancy becomes w — 1/4 — 
7C(3)Z7/8tt 3 -93C(5)C/ 3 /2 9 tt 5 , which is obtained bymak- 
ing the use of energy expansion with respect to U \Lniu\ . 
Thus, the local entanglement near the critical point is, 



E v — 2 — 

In 2 



7C(3)[/ 



n 2 



2tt 3 



Clearly, E v is analytic in the neighborhood of the critical 
point U = 0. 

The variation of the local entanglement caused by 
chemical potential: By adding a chemical potential term 
— rij to the Hamiltonian Eq. l[TJl. the particle num- 
ber of the ground state, hence the filling factor, could be 
tuned. We show the relations between local entanglement 
and the filling factor n for various on-site couplings U in 
Fig. [21 We only need to plot the part for n — N/L < 1 
because the part for n > 1 could simply be obtained by 
the mirror image relation, as easily seen from Eq. JSJ, 
namely 



E v (n)=E v {2-n). 



(11) 



Fig. manifests that the ground state of the half-filled 
band is not maximally entangled as long as U > 0, 
whereas, the maximum of E v lies in between n = 2/3 
and n = 1. There is no double occupancy when U = oo, 
which implies that w = and u + = u~ = N/2L. Hence 
we have an analytical expression of the local entangle- 
ment E v = — (1 — n) log 2 (l — n) — nlog 2 (?i/2) which has 
a maximum at n — 2/3. It is worthwhile to point out that 
at 1/3 filling (i.e., n = 2/3) when U = oo, the ground 
state is a singlet of SU(2|1) Lie supersymmetry algebra 
which possesses the maximal symmetry allowed, while at 
1/2 filling, it is a SU(2) singlet. For U = the ground 
state is invariant under SO (4) rotation at 1 /2 filling. This 
demonstrates that the local entanglement reaches a maxi- 
mum value at the state with maximal symmetry. Accord- 
ingly, the maximum position for < U < oo is expected 
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FIG. 2: The left is the plot of local entanglement E v versus 
the filling factor n for various interaction strength; the right 
is the relation between dE v /dn and U closing to half-filling. 



to lie between n — 2/3 and n = 1, which is confirmed in 

Fig. m 

Except at half-filling where it becomes a Mott- 
insulator, the system is an ideal conductor 

Gl[l3. Con- 

sequently, the local entanglement E v is not smoothly con- 
tinuous at n = 1 for U ^ 0. It is useful to observe the 
derivative of E v with respect to U . 



(10) dE v 
dn 



(log 2 u + - log 2 z) 



n=l~ 



1 dAE 
- + 2 

2 dU 



(12) 



where AE is the gap of charge excitation. Ea. fll|> 
gives rise to dE v /dn\ n=1+ — — dE v / dn\ n=1 - . Obviously, 
there exists a jump in the derivative of E v across the 
point of insulating phase (see Fig. fright) unless U — 0. 

The dilapidation of the local entanglement induced by 
the magnetic field: Since the Zeemann term commutates 
with the original Hamiltonian QJ, the Bethe-ansatz so- 
lution is still applicable. The application of an external 
magnetic field attempts to flip electron spins align along 
the magnetic field. There will be more electrons with 
up spins than that with down spins in the presence of 
magnetic field, consequently, u + > u~ . The spin flip- 
ping will also diminish the density of double occupancy, 
which makes both w and z decrease. As a result, the 
local entanglement is a descending function of the mag- 
netization m z per length. We plot curves of E v versus 
m z for different interaction strength in Fig. [3] The local 
entanglement descends to zero when the magnetization 
is saturated to m z = 1/2. Obviously, the strength of the 
on-site interaction suppresses the local entanglement also 
in the presence of magnetic field. 

^From the above discussions, we find that the local en- 
tanglement manifests distinct features at the point where 
quantum phase transition undergoes. The variation of it 
caused by the chemical potential and its descendence in- 
duced by magnetic field are both not continuous at the 
quantum phase transition points. This is similar to other 
studies, e.g., the one-dimensional XY model in a trans- 
verse magnetic field 0, @i where the derivative of the 
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FIG. 3: The local entanglement E v versus the magnetization 
m z = M z /L for different on-site interactions. 

pairwise concurrence C with respect to the dimension- 
alless coupling constant develops a cusp at the quan- 
tum phase transition point. However, such discontinu- 
ity is not universal. The local entanglement in the one- 
dimensional Hubbard model (for fixed chemical potential 
without external magnetic field) is smoothly continuous 
at the transition point driven by the on-site interaction 
aforementioned. Similarly, we also obtained in our re- 
cent work on the one-dimensional XXZ mode fl2T |. near 
the quantum phase transition point A = 1 where Mott- 
insulator phase transition occurs, the pairwise concur- 
rence C is C = C - Ci(A - l) 2 . No discontinuity occurs 
in these two cases. 

It is indicated recently that two mechanisms 

may bring about quantum phase transitions in one- 
dimensional system. One is caused by the level cross- 
ing of the ground state and the other arises from the 
level crossing of the low-lying excited state where no 
level crossing occurs at the ground state. In the later 
case, the ground-state wavefunction is smoothly continu- 
ous with respect to the variation of parameters that drive 
the quantum phase transition, and hence the entangle- 
ment is also smoothly continuous at the quantum phase 



transition point. For the former case, the level crossing 
of the ground state will clearly cause the entanglement 
to be none smoothly continuous at the transition point. 
Therefore the continuity properties of the first derivative 
of the local entanglement might be an ancillary tool to 
judge the mechanism of quantum phase transition pro- 
posed in [lg| |. 

In summary, we have extensively studied the local en- 
tanglement for the one-dimensional Hubbard model. At 
half filling, we obtained the relationship between the lo- 
cal entanglement and the on-site Coulomb interaction. 
Our result indicated that the local entanglement reaches 
the maximum value at the critical point where Mott- 
insulator transition undergoes, and that the interacting 
strength suppress the local entanglement dramatically. 
Scaling behavior closing to the critical point U = was 
given as manifests that the local entanglement is an an- 
alytical function of U at this critical point. The asymp- 
totic behavior of the local entanglement at the strong 
coupling limit, U — > oo, was also given, Furthermore, we 
calculated the dependence of the local entanglement on 
the filling factor and the magnetization. The variation 
of local entanglement caused by the chemical potential 
shows that the local entanglement reaches maximum at 
filling factor n between 2/3 and 1. For any finite on- 
site Coulomb interaction U, the local entanglement de- 
velops a cusp at n = 1. At the strong coupling limit, the 
1/3 filled band takes the maximum local entanglement, 
suggesting that the ground state with maximal symme- 
try possesses the maximum magnitude of local entangle- 
ment. The local entanglement decreases with increasing 
magnetic field and approaches to zero at saturated mag- 
netization. For all band fillings, with or without mag- 
netic field on, the on-site Coulomb interaction always 
suppresses the local entanglement. 
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